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We study the time evolution of N q two-level atoms (or qubits) interacting with a single 
mode of the quantised radiation field. In the case of two qubits, we show that for a set 
of initial conditions the reduced density matrix of the atomic system approaches that of a 
pure state at ^ , halfway between that start of the collapse and the first mini revival peak, 
where t r is the time of the main revival. The pure state approached is the same for a set of 
initial conditions and is thus termed an 'attractor state'. The set itself is termed the basin of 
attraction and the features are at the center of our attention. Extending to more qubits, we 
find that attractors are a generic feature of the multi qubit Jaynes Cummings model (JCM) 
and we therefore generalise the discovery by Gca-Banacloche for the one qubit case. We give 
the 'basin of attraction' for N q qubits and discuss the implications of the 'attractor' state in 
terms of the dynamics of N q -hody entanglement. We observe both collapse and revival and 
sudden birth/death of entanglement depending on the initial conditions. 



I. INTRODUCTION 

The dynamics of two level quantum systems (also known as qubits) such as spins in a magnetic 
field, Rydberg atoms or superconducting qubits, coupled to a single mode of an electromagnetic 
cavity are of considerable interest in connection with NMR studies of atomic nuclei 111], cavity 

n n 

quantum electrodynamics [2( and the physics of quantum computing [3[. The simplest model that 
captures the salient features of the physics in these fields is the Jaynes-Cummings model (JCM) 
[4], which deals with only one qubit and its generalisation by Tavis and Cummings to the case of 
multiple qubits [fj. 

These models have many interesting features j^], but perhaps the most surprising prediction is 
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the phenomena of 'collapse and revival' of the Rabi oscillations in the qubit system when the field 
is started in a coherent state \a) Q], 

oo n 

\ a)=e -H 2 /2^2 " | n ), a = \a\e- i0 . (1) 



n=0 v 



Here the state |n) is the eigenstate of the photon number operator h = wa with eigenvalue n and 8 
is the initial phase of the radiation field. \a\ determines the size of the field, with n = (a\h\a) = \a\ 2 . 

In short what is found is that the initial Rabi-oscillations of the probability of being in a given 
qubit state decay on a time scale called the collapse time, t c , but then revive after a much longer 
time, t r , due to the quantum 'graininess' of the radiation field [§[ 0, [ijj. Whilst the nature of this 
remarkable phenomenon in the case of one qubit is now well understood [6|, both experimentally 
[3] and theoretically 11]. the multi qubit case 0, [3] has not been explored to the same extent. 



For instance, the occurrence between t c and t r , of what we shall call the 'attractor' state has only 
been fully investigated in the one qubit limit. This intriguing aspect of quantum dynamics was 
discovered and highlighted by Gea-Banacloche in his study of the one qubit case. In this 

paper we report on our investigation of analogous 'attractors' in cases of more than one qubit. 

The organization of this paper is as follows. In Sec. [IIJ the dynamics of the resonant JCM 
(single qubit coupled to a cavity system) is presented. This sets up our notation, analytic tools 
and methods of data presentation. In Sec. IIIII we extend the model to two qubits and show that 
the idea of an 'attractor' state still exists. We investigate the dynamics in depth and quantify the 
set of initial states that approach the 'attractor', giving the 'basin of attraction'. We show that the 
atomic system exhibits interesting entanglement properties, in particular that the entanglement 
can vanish and then return at a later time. Furthermore, the manner in which this occurs depends 
on whether we are in or outside this basin of attraction. In Sec. IIVI we develop these ideas further 
from those of the two-qubit case and find a set of initial states (the basin of attraction) from which 
the system always goes to the 'attractor' state for an arbitrary number of qubits. In the limit of 
large qubit number we rewrite our qubit equations in terms of spin coherent states and demonstrate 
some intriguing phenomena associated with the iVg-qubit JCM, especially the collapse and revival 
of A^-body entanglement. Whilst our method of analysis, and some of our results, are similar to 



those of Chumakov et al. I15H and Meunier et al. 



la ], our focus is on the quantum information 



aspects of the problem and in particular the dynamics of entanglement. 



3 



II. FEATURES OF THE JAYNES-CUMMINGS MODEL 



lzi, HQ. 



We first review the JCM and appearance of 'attractor' states in the one qubit case 
The JCM model consists of a single qubit coupled to a single mode of a quantised electromagnetic 
field. The two possible states of the qubit are the ground state \g), with energy e g , or its excited 
state |e), with energy e e . We consider the following Hamiltonian: 

m - ■ - • ~ + -^ (2) 

(3) 



H = hwa)a H a z + h\ (aa + + a) a 

2 V 

o z = |e)(e| - \g)(g\, a + = \e)(g\, a~ = \g){e\ 



where a) and a are the creation and annihilation operators of photons with frequency u, HQ, = £ e —£g 
and A the cavity-qubit coupling constant. The model is only valid close to resonance and we shall 
only consider the resonant case of uj = O. 
We have an initial state: 



|*i(t = 0)) = = (C g \g) + C e \e)) ^ C n 



n 



(4) 



n=0 



1 2 1 2 

where the state of the qubit is normalised so \C e \ +\C g \ =1 and the field is started in a coherent 
state. In the interaction picture (i.e. in a rotating reference frame) we obtain a solution to the 
Hamiltonian: 



= ^2 [{CeC n cos (\Vn + It) - iC g C n+ i sin (\Vn + It)) \e, 



n) 



n=0 



(C g C n+ i cos (AVn + It) - iC e C n sin (AVn + It)) \g, n + 1)] + C g C \g, 0). (5) 

In spite of its relative simplicity, Eq. ([5]) exhibits a wide range of interesting phenomena. A 
much studied example is the collapse and revival of Rabi oscillations. If the qubit is started in the 
ground state then initially the probability of the qubit being in the ground state oscillates from 
one to zero at the Rabi frequency. These oscillations then decay on a time scale called the collapse 
time, t c = j, but then revive after a much longer time: 



A • 

When n > At we can approximate in Eq. ([5]) by ^i(i)): 



(6) 



+ ±{t) D,i(t)) $ + i(t)) + /Li(t) D_i(t)) <i?_i(m, 



(7) 
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where /3 ± i(t) is a normalisation factor, D±i(t)) is a state of the qubits and <&_j_i (t) \ is a state 



of the field. 



e w C e =f Q 



±1 



e Te 



^fi7rt/t r I 



(8) 

(9) 
(10) 



With the solution written in this form it is clear there are two distinct parts to the wavefunction. 
In fact the field just consists of two different coherent states that are out of phase by Xt/^/E. In a 
penetrating study of this model Gea-Banacloche [lil noted that at a time %, in the very large n 
approximation, the qubit disentangles itself completely from the field and the wavefunction factors 
into a product state of the radiation field and the qubit. 

To make this more apparent we rewrite Eq. ([< 



Eq. (HH at the time t = f. 



± 1 c ±i-Kn/2 

V2 



(I) 



i 

\±ia). 



With the wavefunction written in this way it can be seen that the two states 
D_i (tj-)/ are in fact the same at the time We label this state 



D, 



Yattl 



ktatt) 



1 

71 



e- ie \e)±i\g) 



(11) 
(12) 
(13) 

and 
(14) 



This state is independent of the initial coefficients C e and C g and only depends on the phase 6 
of the initial coherent state, Eq. ([1]). Throughout this paper we follow Phoenix and Knight [20I 
and call this state an 'attractor' state. Note that we use the term 'attractor' to refer to (V'latt) 
even though in the theory of non-linear differential equations the term attractor is reserved for 
solutions which, when reached, after evolving from an arbitrary initial condition, stop changing 
with time. As a consequence the qubit state can be factorised out of the wavefunction meaning 
the qubit and the field are in a product state: 



ie 



iim/2 [ (J 



a 



ia) 



ie 



-i-wn/2 



e l9 C e + a 



-ia, 



(15) 



It is clear that all the information of the initial state about the qubit is stored in the radiation 
field which happens to be in a "Schrddinger cat" state. At a later time t = ^k 1 , the qubit again 
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disentangles itself completely from the field and goes to a pure state. The qubit state at this 
time is orthogonal to \^t a tt) anc ^ we l aDe l it IV'fatt)- This remarkable behaviour is most unlike 
the consequences of a linear Schrodinger equation satisfied by a qubit without the cavity field or 
with a classical field and, as was stressed by Gea-Banacloche, is a natural route to quantum state 
preparation. All initial qubit states tend to the same 'attractor' state (which is determined by 
the initial phase of the field state) and the quantum information about the initial qubit state has 
effectively been swapped out and encoded into the field at the 'attractor times'. 

Before we move on to describe the results of numerical simulations we consider the field, Eq. 

$, become 



(jlOp . at two revival times t = t r and t = 2t r . At both these times the field states 
equal, with = < &„i(t r )^) = |— a) and $ + i(2f r )\ = <l>_i(2t r )^ = |a). As a consequence 

the system is once again in a product state as the field can be factorised out of the wavefunction 
for the whole state. The information about the initial state of the qubits has been transferred from 
the field back to the qubits. 

By considering the average number of photons in the field, n, to be large many interesting 
observations have been made of the JCM. It was shown that the system repeatedly evolves from 
a product state to an entangled state and back again. Although the analytical equations for large 
values of n clearly show this behaviour, we now consider the exact solution \^i(t)) for modest n 
(n = 50), shown in Eq. ([5]) and comment on the results of numerical simulations for the qubit and 
the photon sector separately. 

We consider the state of the qubits after tracing out the field state. The probability of the 
qubits being in the state \g) is, 

oo 

P lg (t) = (g\p q (t)\g) = Y,\(g,n\* 1 (t))\ 2 (16) 

n=0 

where p q (t) = Trj (\Px(i)|) is the reduced density matrix at a certain time for the qubits 

when the field has been traced out. Throughout this paper all instances of the letter q refer to 
the qubit system and / refers to the field. The results are depicted in Fig. [T] for the initial state 
\^i(t = 0)) = |fl r )|a) and the well known phenomenon of collapse and revival can be seen. For 
numerical reasons the sum over n is truncated at n max = 200 as this is adequate for the average 
number of photons in the mode of n = 50. 

An entropy measurement of the system can be made at any point in time. Entropy is a well 
defined quantity that measures the disorder of a system and the purity of a quantum state. In cases 
where the time-dependent Schrodinger equation (as opposed to some master equation for reduced 
dynamics) determines the time evolution of a system, the value of the entropy remains constant. 
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FIG. 1: (color online) Time evolution for a system with one qubit when the qubit starts in the initial state 
\g) and the value of n = 50. (a) the entropy of the qubit. (b) the probability of being in the qubit's initial 
state \g). (c) the probability of being in the state \ipatt)- At 4p the probability of being in the 'attractor' 
state goes to one while the entropy goes to zero. 

The complete systems (qubits plus field) we study are always started in a pure state and follow 
Schrodinger dynamics; therefore the entropy of the complete systems always remains zero as we 
have a pure state. In these systems the individual entropies, of the qubits or the field, are of more 
interest. We calculate the partial entropy from the reduced density matrix p q using the rescaled 



von Neumann entropy 



21 



228 



S q (t) = -Tr(p q (t)log 2Pq (t))/N q , 



(17) 



which is rescaled such that the entropy of the qubits ranges from zero, for a pure state, to one, 
for a maximally mixed state. For the bipartite case (the one qubit N q = 1 JCM) this reduces to 



the standard von Neumann entropy. As was discovered by Araki and Lieb [23], if the system is 
started in a pure state then S q (t)=Sf(t) for all time, where Sf(t) is the entropy of the field system 
calculated from reduced density matrix pf(t) = Tr q Therefore for the rest of the 

paper we shall refer to the value of S„(t) only when we mention entropy. 

In quantum information these partial entropies provide a measure of the entanglement present 



between the two systems, known as the entropy of entanglement [2j, |25(. If the entropy is zero 
then the state being measured is pure and there is no entanglement present. If the value of the 
entropy is non zero then there is some entanglement present between the two parts of the system 
and the reduced density matrix is mixed. Returning to the JCM we note that the 'attractor' state, 
Eq. iTH]) . is a pure state and hence the value S q (^-) will be zero. So the qubit has disentangled 
itself from the field and there is no entanglement present at this time. The reduced entropy is a 
powerful tool in understanding entanglement in the JCM and it is worthwhile recalling an example 
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of its power to elucidate the quantum dynamics [20 . |26 ] . 



We have numerically evaluated S q for qubits started in the initial state \^\{t = 0)) = 
using the exact solution for and the results are shown in Fig. [l]by the curve labelled (a). 

We also plot the probability Pi a tt(t) = (4>tatt\p<i(t)\' l l ; iatt) that the qubit is in the 'attractor' state 
\ipi att ), as defined in Eq. (|14H . shown by the curve labelled (c). 

Evidently at the time y the reduced entropy goes very close to zero and Piatti^f) goes very 
close to one. As mentioned above, this means there is very little field-qubit entanglement present 
and the qubit is essentially in the 'attractor' state, IV'iatt)' as defined in Eq. (|14p . At ^ the 
probability of being in this 'attractor' state approaches zero, as the qubit is in the state \ipi att ), 
which is orthogonal to | iptatt ) • These results all agree with the discussions of Eq. (fT5|) . At future 
'attractor' times the system alternates between |V , i~ a tt) ana - l^iatt)- 

A big difference between numerical results and the large n solution is noticed at the time t = t r . 
Earlier we showed that the solution in the large n approximation will return to a pure state at this 
time, but in Fig. Q]the entropy has a large value of 0.7, wedged between two maxima (indicating 
large entanglement between qubit and field). This is because the width of this dip is on a time 
scale similar to the collapse time, much narrower than the entropy dip at y. As the value of h is 
increased in numerical simulations, the value of the entropy at t = t r gets closer to zero, but we 
see that even for fairly large coherent states, the large n approximation can deviate significantly 
from the full solution. 

So far we have reviewed the time evolution of the reduced qubit density matrix in the one-qubit 
and one mode system. Before we move on to address the two-qubit case we calculate the Q function 
[27I to investigate the radiation field, 

Q(a,t) = (a\p f (t)\a). (18) 

This is shown in Fig. [2] in the rotating reference frame for the exact solution as contour plots. 
We give an example of a 3d plot of Fig. [2](b) in Fig. El but for the rest of the paper we will use 
the 2D contour plots. 

The Q function shows a set of two wavepackets of the cavity field, each of which can be thought 
of as a Gea-Banacloche state [ljj] in the large n solution. As a function of time the two 'blobs' 
move around the complex plane and follow the circular path of radius \[Tl. These 'blobs' start to 
smear out and become distorted from the original coherent state. Although the wavepackets all 
begin in the same place, they evolve with different frequencies, so the states begin to separate. 
After a period of time, depending on the differences between the frequencies, the different 'blobs' 
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FIG. 2: (color online) A contour plot of the Q function at six different times when the qubit is in the initial 
state \g) and n = 50. The atomic dipole states are represented by arrows, (a) the time t = 0. The system is 
started in a coherent state which is shown by a circle of uncertainty in phase space, (b) a time a little after 
that of figure (a), t = The 'blobs' are well separated and are moving around a circle of fixed radius, (c) 
the time t = %. This corresponds to the first dip in entropy, (d) just before the time t = t r . This is just 
before the revival time and the two 'blobs' will overlap, (e) the time The second dip in the entropy, (f) 
the time t = 2t r . Both the wavepackets have returned to their original position. 
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FIG. 3: (color online) The Q function in 3D for the time t = ^ corresponding to diagram (b) in Fig. @. 
We notice two distinct peaks of equal height that are very similar to coherent states. 



are separated by more than their diameter and can be easily distinguished. Eventually these two 
wavepackets will again overlap and this occurs at the revival time, t = t r . When these 'blobs' 
have some overlap in phase space there will be Rabi oscillations in the qubit probability, giving the 
periodic sequence of revivals seen in Fig. [TJ However, as can be seen in Fig. EJd), at the revival 
time, t r , the 'blobs' have undergone some distortion and so the wave packets are no longer simple 
coherent states. This distortion elongates the wave packets, and so they overlap for a longer time 
than in the ideal coherent state case. We therefore understand the lack of a full amplitude revival, 
Pig(U) = 1, as due to the spreading, and hence distinguishability, of the field states. When the 
value of n gets very large there is less distortion and the Gea-Banacloche states stay exact coherent 
states throughout their evolution. 

Also included in Fig. [2] is an arrow for each 'blob' which represents the average atomic po- 



larisation, di i(t) 



a 



D,i(t)\ in the equatorial plane of the Bloch sphere. Fig. Wi c ) 



shows the two dipole arrows are in line and pointing in the same direction at the time y because 



D-l{^r)) = Intuit)- The arr ows in Fig. [2]|e) are again both the same, but are 



pointing in a different direction as 



D 



'3t r 



D 



'3t r 



at I 



>. The field states are mov- 



ing in a different direction in Fig. [2je) to that of Fig. EKc) and their identity has swapped. These 
arrows can be used to see if the two qubit states are the same and the direction of the dipole arrows 
highlights we are in a different 'attractor' state. 

Inspired by the above behaviour, in this paper we extend the analysis and investigate new 



10 



problems. In particular we address the question: "Are there similar 'attractor' states in multi 
qubit systems and what are the implications of their existence?" 



III. TWO QUBIT JAYNES-CUMMINGS MODEL 

The N q qubit Hamiltonian, which in this section we treat for N q = 2 so i = 1, 2, is a generali- 
sation of Eq. ([2|). It is called the Tavis-Cummings Hamiltonian 



H = hwa ] a + - ^2 fli&i +h^2\i (aaf + a) dr J . (19) 

i=l i=l 

For brevity we only consider the solution for the case Ai = A2 = A and uo = VL\ = This model 



has been studied in detail by Chumakov et al. 



and the large n solution has been found Q] 



with similar results to our analysis. Adding another qubit to the system yields a more complex 
system which can exhibit interesting entanglement properties as not only can the two qubits be 
entangled with the field, but they can also be entangled between themselves. 

We have investigated the evolution of the system over time, with analytical solutions in the 
large n approximation [29J and the numerically evaluated exact solution. To summarise these we 



start with the initial state: 



|*2(* = 0)) = \ifo)\a) = (C ee \ee) + C eg \eg) + C ge \ge) + C gg \gg)) \a) 



(20) 



where the state is normalised so |C ee | + \ C eg \ + \ C ge \ + \ C gg \ = 1. Following the same procedure 
as in the one qubit case we solve the Hamiltonian and approximate for large n. The time evolution 
of this state takes the following form: 



*2(*)) = Yl k (t)\D k (t)) ® |* fc (t)>, 



k=-l 



where k = —1, or 1 and: 

±2i7r-^(n+2) 



0±l(t) 

00 

|A)(*)> 



(e 2 * e C ee =p e w (C eg + C ge ) + C gg ), 



1 
\D 



f>2idr< _ n \ l \\n _ c, I 2 



1 



C gg ) {e-™\ee) - \gg)) + (C eg - C ge )(\eg) - \ge)) 

2 ^|e 2 *^C ee C gg \ |Ce<; I 
T2i7Tt/t rp -i0 (\ pn \ , \ np \\ , p T4i7Tt/i 



e-^(|e 9 ) + | 5 e)) + e 
D ± i (2i)) ® D ±i (2i)Y 



l<?<7>], 



!**(*)> 



e 2ik-Kt/t ra 



(21) 

(22) 
(23) 

(24) 

(25) 
(26) 
(27) 
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Note that the state \Do(t)) = \Dq) is time-independent and is thus a maximally entangled qubit 



state for all times if C, 



'■9 



a 



y<-' 



29l | . The field states are called Gea-Banacloche states and once 



again have a specific qubit state associated to each of them. 

By making the observation shown in Eq. (|26p it is easy to see that the two qubit states \D±i(t)} 
are direct products of the one qubit states D ± i(t)\ up to a mapping t — ► 2t. This can be used 



to predict the behaviour of \D±i(t)}. As stated in Sec. [UJ D, i(%) 



att/ 



conclude that \D±\{t)) will go to a direct product of V'latt) a * * 



l • 



D±\ ft) 



matt) ® I V'latt) 
I ^att) 



therefore we can 

(28) 
(29) 
(30) 



So we have found a term that is similar to the one qubit 'attractor' state, V'latt)) w hich we 
shall label the two qubit 'attractor' state, (V'Jatt)- 

1 



2 att) 



e 2i »\ee)±ie™(\eg) + \ge))-\gg) 



(31) 



In the one qubit case the qubit system went to a second 'attractor' state | V'latt) a ^ the time ^ 
and means |-D±i(^)) = |^att) 8) |Vw) = \^2att)- 



Basin of Attraction 



The two qubit JCM is more complicated than the original one qubit case and shows some very 
interesting properties. Firstly we rewrite the wavefunction for the whole system as the time t = ^ 
to understand the two qubit case in more detail. 

1 



" o I V'latt / 



2ie 



C e e "I" C, 



e iS {peg + Cg e ) ( e 



Jnn/2 



99 



ia) 



e i7rn/2 \ia) + e- inn/2 \-ia) 



+ Po\D )\a) 



(32) 



With the wavefunction written in this form we can see which initial conditions result in the 
entropy going to zero. Indeed when (3q = the qubit and field parts are each in a product state, 
so there is no qubit-resonator entanglement present. There is also no qubit-qubit entanglement as 
these are in the unentangled 'attractor' state and this means there is absolutely no entanglement 
in the system at t = Interestingly, the information on the initial state of the qubits is again 
stored in the radiation field and the qubits are in a spin coherent 'attractor' state (see appendix 

ED. 
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The initial conditions of the qubit that results in (3q = are: e l9 C ee = e td C qq = a and 



C eg = C ge = \j \ — \a\ 2 . So the initial state of the qubits will be of the form: 



|V> 2 > = a(e- i9 \ee) + e i6 \gg)) + \a\\\eg) + \ge)) (33) 

where 9 is the initial phase of the radiation field and < \a\ < These states lie in the 

symmetric subspace and have some interesting properties that will be discussed in section IIII Bl 
While the qubit state is within this basin of attraction its state is parameterized by the single 
complex variable a and the wavefunction only contains two field states, |<3>-i-i(i)). At the time ^ 
these two field states are macroscopically distinct coherent states: 



*2(|)) = ~e W \^ 



2att/ 




inn/2 [ / 1 I |2 1 I • \ , -inn/2 I . /I i |2 . 

I a + \l - — \a\ I \—ia 



(34) 



We now see that the initial state of the qubits alters the relative proportions of \a) with respect 
to |— a). This is also seen in the one qubit case and is the mechanism whereby the quantum 
information of the initial qubit state is encoded in the field. State preparation of a Schrodinger 
cat state utilizing this phenomenon has been discussed [14I I30I], With the two qubit case the same 
protocol can be used, but is not so useful as the initial qubit state has to be within the basin of 
attraction in contrast to the one qubit case where all initial qubit states lead to a Schrodinger cat 
state at the 'attractor' times. 

There is also a connection between the field states \&±i(t)} in the two qubit case and the 
field states $>±i(t)^ from the one qubit case, |$±i(t)) = <3?_|_ 1 (2i)^ . Again there is a mapping 
t — > 2t. In Sec. [TT] it was shown that the field and qubit become unentangled at t r and 2t r as 
$ + i(t r ))= $_i(t r ))and $_|_i(2t r ))= < &_i(2i r )\. Using this information from the one qubit 



2 

case we can predict that |$ + i(y)) = |<3?_i(y)) and |<I>-|_i(t r )) = |<I>_i(t r )) which is indeed the 



case. As a consequence the field part of the wavefunction can be factorised out and there is no 
entanglement present between the qubits and the field at the times y and t r . 

We now consider the numerical results to test the analytical predictions and see how the system 
behaves when n = 50, shown in FiglU We consider the probability of the qubits being in the state 
\gg) and the probability P2att(t) = ( 1 p2att\Pg(^)\' > p2att)^ ^ na ^ l^tatt) occurs a * a particular time for 
two sets of initial conditions, one satisfying (3q = (Figl^b)) and one not (Fig[Ha)). Evidently the 



pattern of collapse and revival is more complicated than in the one qubit case 



31[, but nevertheless 



we still see characteristic dips in the entropy. The value of the entropy at is not the same in 
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the two cases shown in Fig. U] and in fact they show two completely different behaviours. In Fig. 
IHa) we notice that the dip in the entropy goes to a value of 0.35, very different to the value of 
the entropy in Fig. [U but in Fig IHb) we do indeed find that the entropy is very small at -f . The 
initial state used for Fig. [U^a) is not in the basin of attraction, whereas the initial state used for 
Fig |T{b) is. 

A second dip in the entropy can be seen in Fig. 0] at the time t = ^ which, in the large n 
approximation, goes to zero. The probability of being in the 'attractor' state \1p2att) g° es to zero 
at this time. The field has disentangled itself from the field again, so at this point the qubits are 
again in a pure state, the orthogonal 'attractor' state j^fatt)- 

It can be seen from the analytic solution that, in the large n approximation, the field and qubits 
also disentangle themselves at the time t r , for all initial conditions of the qubits, and also at y 
if the initial qubit state is in the basin of attraction. This is not seen in the diagrams in Fig. U] 
and instead the value to the entropies at these times are fairly high which means the numerically 
evaluated exact system is far from the analytic results in large n approximation. This is similar to 
the one qubit case and can be explained by considering the dynamics of the field. 

We can once again plot the Q function for the numerical solution to pictorially display what is 
occurring in the cavity field. The values of Q{a, t) in the complex a plane at fixed times are shown 
in Fig. [5] when the qubits are initially in the state \gg), which is outside the basin of attraction. 
This time the Q function shows a set of three wavepackets of the cavity field (compared to the 
two for the one-qubit case, Fig. [2]). As a function of time the three 'blobs' separate, but only two 
move around the complex plane and follow the circular path of radius Again the arrows show 
the dipole moment of the qubit state associated with each wavepacket, dk{t) = {Dk{t)\ J2\Dk(t)). 
The 'blob' that does not move over time has a dipole moment of zero, so only the two moving 
Gea-Banacloche states have an arrow. 

For the one-qubit case the revival time occurred when the two wavepackets overlapped. In this 
case we have two 'blobs' that move around phase space and overlap at the time t = ^, leaving the 
stationary wavepacket at the initial position. This causes a rather weak revival peak as only two 
of the three wavepackets are in phase. At a later time, Fig. all three wavepackets converge at 
the starting position and there is a main revival peak, t = t r . If the initial state of the qubits was 
in the basin of attraction then (3q ~ and there would only be two wavepackets. The Q function 
pictures would then look similar to those in Fig. [2] and the system would behave like the one-qubit 
case. Nevertheless we note that the radiation field returns to the start position after one revival 
time, t r compared to the one-qubit case where it returns to the same place after 2t r . 
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FIG. 4: (color online) Time evolution for a system with two qubits when n = 50 for two different initial 
states of the qubits. (Top)The initial state of the qubits is \gg) and it can be seen that the two qubit 
'attractor' state is not reached at ^. (Bottom) The initial state of the qubits is ^(|ee) + \gg)) (i.e. within 
the basin of attraction) and it shows that the two qubit 'attractor' state is reached at 

In both diagrams (a) the entropy of the qubits. (b) the probability of the two qubit state \gg). (c) the 
probability of being in the two qubit 'attractor' state \1p2att) wrien the initial phase of the radiation field is 
6» = 0. 

In summary when there are two qubits interacting with a single cavity mode the wavefunction 
consists of a part that behaves like the one qubit case, \D±i(t)) and an extra part \Dq). The 
'attractor' state, \ipf att ), occurs at the times when the two 'blobs' in phase space are macroscop- 
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FIG. 5: (color online) The Q function at six different times when the qubits arc initially in the state \gg). 
Once again the atomic dipole states are represented as arrows and n = 50. (a) the time t = 0. The system 
is started in a coherent state which is shown by a circle of uncertainty in phase space, (b) a time a little 
after that of figure (a), t = ^. (c) the time t—if. (d) the time t = (e) the time (f) the time 
t = t r . Both the wavepackets have returned to their original position. 
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ically different and when (3q = 0. Moving from the one qubit case to the two qubit case involves 
a mapping of t — > 2t and a similar mapping can be used to understand the equations for more 
qubits. 



B. Collapse and Revival of Entanglement 

An interesting new feature of the two-qubit case, as opposed to the one-qubit case, is that there 
can be entanglement present between the two qubits themselves, not just between the qubits and 
the field. So far we have used the entropy to measure the amount of entanglement between the 
field and the qubit system as the whole system is in a pure state. Unfortunately this measure can 
not be used to calculate the entanglement between the two qubits as they are sometimes in a mixed 
state themselves. To find the amount of entanglement between the two qubits we use the mixed 
state tangle r, defined as [32. l33f| : 



r = [max{Ai - A 2 - A 3 - A 4 , 0}] 2 . (35) 

The X's are the square root of the eigenvalues, in decreasing order, of the matrix 
Pq(cr y <8> cr y ) p* q (a y ® cr y ), where a y = i(a~ — a + ) (see equation Eq. [3]) and p* represents the 
complex conjugate of p q in the |e), \g) basis. A tangle of unity indicates that there is maximum 
entanglement between the two qubits. If the tangle takes the value of zero, then there is no en- 
tanglement present between the two qubits. To clarify, the entropy will be used to measure the 
amount of entanglement between the field and the qubit system whereas the tangle will be used to 
measure the entanglement between the qubits. 

In this subsection we focus on the initial states inside the basin of attraction. The value of 
the tangle for the states given by Eq. ([33]) can be seen in Fig. [U This diagram shows only two 
points where r = 0, for a = ±5, which demonstrates that there are only two product states in 
the 'basin of attraction'. In contrast, there are many points in which r = 1, for a = ^7= where d> 

v2 

is an arbitrary phase, or a = ir where r is a real number. This means that all possible levels of 
entanglement are represented in the basin of attraction, Eq. (|33j) . 

Whilst almost all states in the basin of attraction given in Eq. (J33]) describe entangled states, 
these all evolve into \iptatt) at t = ^ and \1p2att) at i = ^l 21 , at which times they are pure and 
unentangled. So if the system is started off in one of these entangled states defined by Eq. (133 j) . 



at the time ^- the system goes to a state which has zero entanglement. Therefore we have gone 
from having entanglement in the system to having absolutely no entanglement, either between 
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FIG. 6: (color online) The value of the tangle for the states in the basin of attraction, for different values of 
a. We notice that the plot shows only two points where the tangle is zero. 



the qubits and the resonator or between the qubits themselves. At a later time this entanglement 
returns to the system and can be said to 'revive'. The tangle returns to its initial value at y and t r , 
regardless of the phase of the coherent state. This can be seen when we consider the wavefunction 
for the system at this time: 



1 



-a) 



e l7Tn a 



e^ e \ee} + \eg} + \ge}+e te \gg) 



+ e~ 



a + M- 



e-^\ee)-\eg)-\ge)+e%g) 



(36) 



As stated in Sec. IIII Al this state shows no entanglement between the qubits and the field and 
therefore the entropy of the reduced qubit density matrix will be zero. The value of the tangle 
for the qubit state at t = y is found to be exactly the same as the tangle for the initial state 
of the qubits, therefore r(0) = t( ^ ). We have called this phenomenon 'Collapse and Revival of 



Entanglement' 



341 ] . Analogous to the single qubit case, the quantum information in the initial 



state of the qubits (for states within the basin of attraction) is swapped out and encoded into the 
field state at the 'attractor' times. Calculations demonstrate that the value of the tangle remains 
near zero for long periods between revival. Therefore we are dealing with a phenomenon which 



has been described as the 'death of entanglement' 



centre of much current interest [37j] . Qing et al. 



by Yu and Eberley 



35. 



361 ] which forms the 



have found a similar collapse and revival for 



the same model we have studied here but for very different initial conditions. Even just within 
the system we consider here, detailed studies reveal interesting links between the dynamics of the 
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disappearance (and reappearance) of entanglement and the basin of attraction, as we discuss in 
this and the subsequent subsection. 

So if the qubit system starts in the basin of attraction and r / 0, over time the initial entan- 
glement in the qubits is exchanged for entanglement between the qubits and the field, and then 
for a superposition of field states. At the time of the 'attractor' state, there is absolutely no 
entanglement present in the system. Then the process starts to reverse. Firstly there is again 
entanglement between the qubits and the field and at the two qubit revival time, r returns to 
its initial value and only qubit entanglement is present. 

The actual extent of the entanglement revival for the exact solution can be seen in Fig. [7J (a) 
where we show the tangle for different initial states within the basin of attraction. On first glance 
we see that the amount of tangle present at t = y is nowhere near the initial amount of tangle 
present at t = 0, which is far from the prediction in the large n approximation. The value of the 
entropy at y is around 0.5 for the case when a = -^=, rather than the predicted value of zero, 
Fig. [7J (b). In order for there to be maximum possible entanglement in the system there has to be 



the least amount of entropy, so there is a trade off in tangle for entropy 

Another very interesting feature illustrated in Fig. [7J (a) is the manner in which the entanglement 
vanishes. For all states within the basin of attraction, although the initial entanglement decays 
rapidly (before reviving at y and later times) it does so with a Gaussian envelope (as per the Rabi 
oscillation collapse) . It therefore goes smoothly to zero — this is reflected by the further observation 
that the eigenvalue expression in Eq. (|35p never goes negative, so the max operation never needs 
to be implemented for the tangle. This can be explained by considering the rank of the matrix 
used to calculate the tangle, p q (a y ® a v ) p* (a v <8> cr y ). The rank of a matrix is equal to the number 
of non-zero eigenvalues. The matrix used to calculate the tangle for this system has a maximum 
of rank= 2 when the qubits are started in the basin of attraction, so A3 = A4 = 0. Due to the 
definition of the tangle, Ai > A2 and Eq. ([35]) becomes: 



T = (Ai - A 2 ) 2 > 0. (37) 

This characteristic 'Collapse and Revival' behaviour is to be contrasted with behaviour where 
entanglement vanishes or appears with finite gradient and where the max operation is needed in 



Eq. ([35]) because the eigenvalue expression does go negative. 



b distinguish these behaviours, we 



use the terminology 'sudden death/birth' of entanglement |35|] for the latter. This phenomenon 
does not occur for initial qubit states within the basin of attraction (which all show collapse and 
revival of entanglement) but can for those outside, as will be demonstrated in the next subsection. 
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FIG. 7: (color online) (Top) 'Collapse and Revival of Entanglement' when the qubits are started in the 
maximally entangled state (|ee) + \gg))/y/2 where n = 50. (a) the entropy of the qubits. (b) the probability 
of the two qubit state \gg). (c) the tangle for the system. 

(Bottom) 'Collapse and Revival of Entanglement' for many different initial qubit states in the basin of 
attraction with n = 50. These are shown in many different colours depending on the initial entanglement 
measure. It shows that the entanglement in the qubits is initially lost then returns at the first revival peak, 
t = % . As the initial state is in the basin of attraction we know the whole system has zero entanglement 
at t = if. So there is even no entanglement between the field and the qubits. Therefore we know that over 
time the amount of entanglement varies a great deal. 
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C. Outside the Basin of Attraction 

The previous subsections discussed interesting features of the entanglement dynamics for the 
two-qubit JCM when considering initial conditions inside the basin of attraction. In this subsection 
we highlight interesting features that emerge when considering all other cases — when the qubit 



Hubert space is not restricted and, for example, the terms (b\eg) + y 1 — \b\ 2 \ge))/y/2, where b is a 
complex number, are present in the initial state, as well as unequal superpositions of \gg) and |ee). 
To date, in other works it has been found that the JCM for two qubits can exhibit an interesting 
phenomenon called 'Sudden Death of Entanglement', for certain initial conditions of the qubits. 
Some investigations start with the qubits and the field in a mixed state and some others do not 



38] 



In our work we consider the case when the radiation field is started in a coherent state. Qubits 
states inside the basin of attraction demonstrate collapse and revival of entanglement, but not 
sudden death or birth. However, if we venture outside the basin of attraction and include other 
initial conditions of the two qubit system, we find that there is indeed sudden death of entanglement. 
This means that we do not restrict the condition /3q = in Eq. (|23[) . 

Now, without this restriction, the qubit system will not go towards the 'attractor' state at ^ 
and ^ in the large n approximation as (3q ^ 0. Furthermore the entropy will be non zero at 
these times as the qubit state can not be factorised out of the wavefunction. However, the system 
will still have zero entropy at the 'full' revival time t r as |$_i(i r )) = \&o(t r )) = |<!> + i(i r )) = \a) 
and the field state can be factorised out of the wavefunction. Like before r(0) = r(t r ) so the 
amount of entanglement between the qubits is the same at the main revival time and t = 0. For 
n = 50, this revival will not be complete and so the tangle will not return to its initial value. 
Fig. [8] shows examples of the qubit entanglement dynamics for initial states which are outside the 
basin of attraction for the numerically evaluated exact solution. In this figure a new measure of 
qubit entanglement has been used called the concurrence £ = ypr and has been used as a direct 



comparison to previous work 



36, 



371, 



Remarkably, Fig. [8] shows that a field initially in a coherent state can give sudden death and 
birth of entanglement, with initial qubit states of various forms lying outside the basin of attraction. 
A case where the initial conditions are in the basin of attraction is shown for reference. Clearly 
peaks in the entanglement measure at both ^ and ^ arise (with death/birth at finite gradient) — 
times when there would be no entanglement if the initial qubit state was in the basin of attraction. 
There are still peaks at the times % and t r , but these have a lower maximum and also now exhibit 
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FIG. 8: (color online) (Top) Values of the concurrence over time for several different initial qubit states. As 



in previous figures, n = 50. (a) The initial state of the qubits is -^=|ee) — y yjj \gg). At times the concurrence 
is greater than the initial value, (b) The initial state of the qubits is (\eg) +i\ge)), and unlike initial 
conditions within the basin of attraction, leads to a concurrence value of zero at 4^. (c) The initial state 
of the qubits is — |=|ee) + u~^j\gg) an d shows 'Sudden Death of Entanglement', (d) The initial state of the 
qubits is -j= (\eg) +e 1 f|ge)^, a maximally entangled state that exhibits 'Sudden Death of Entanglement', 
(e) The initial state of the qubits is i (|ee) + \gg}), a maximally entangled state that lies in the basin of 
attraction that is used as a comparison. 

(Bottom) 'Sudden Death of Entanglement' when the qubits are started in the state -^=|ee) + y|§| <?<?)) case 
(c) in the top diagram. As before (a) the entropy of the qubits. (b) the probability of the qubit being in 
the state \gg). (c) the mixed state concurrence of the qubit system. 
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birth/death behaviour rather than Gaussian collapse and revival. One initial condition that yields 
'Sudden Death of Entanglement' has been shown in Fig. [8] alongside an entropy measurement. 

To recap, we find that our system can display both a sudden birth/death of entanglement as well 
as a Gaussian collapse and revival of entanglement, depending on initial conditions. Furthermore, 
we understand that, in this system, the sudden death occurs when Ai — A2 — A3 — A4 becomes 
negative and so, according due to the max operation used in Eq. [35] the tangle goes to zero with a 
finite gradient. In contrast, within the basin of attraction, there are only two non-zero eigenvalues 
Ai > A2, so the tangle can only display a smooth collapse. 

These additional peaks in the qubit entanglement seen in Fig. [Hare due to the inclusion of \Dq) 
and |3>o(£)) i n the wavefunction. From the analysis in Sec. IHII we showed that |-D±i(^)) = \1P2att)' 
the unentangled two qubit 'attractor' state, so all entanglement present at t = is due to |.Do)- If 
I .Do) is a product state then there will be no entanglement in the qubit system. However if \Dq) is 
entangled it still does not mean there will be some qubit entanglement as a mixture of entangled 
and product states need not be entangled. For example curve (b) in the top diagram of Fig. [8] 
shows there is zero concurrence at ^, even though \Dq) contains some entanglement. 



IV. N q QUBIT JAYNES-CUMMINGS MODEL 



Having found a sense in which the one qubit 'attractor' state has analogues for two qubits, it is 
natural to inquire whether similar generalisation holds for the case of more than two qubits. We 
use the Hamiltonian shown in Eq. (|19p for N q qubits. We have found that there are indeed cases 
where the system goes to a state similar to the one qubit 'attractor' state and have a general form 
for the basin of attraction for N q qubits. 

In the large n approximation the state of the system takes the following form (l^ ]: 



N„/2 

Pk(t)\D k (t))®\* k (t)). 

k=-N q /2 



(38) 



|-Dfc(t)) has a complicated form and occurs in the paper by Meunier et al. [if]]. For this work we 
will be most concerned with the states that appear in the expression for the basin of attraction. 
These take the following simple form: 



^ 2 



(t) 



D ±h (N q t) 
D ± UN g t) 



(39) 
(40) 
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Note that: 

|$ fc (t)) = ^ kN ^ a \ (41) 

and we can use this equation to understand the simple dynamics of the N q solution. The field 
initially starts in a coherent state that then splits into N q + 1 different 'blobs' and a revival 
peak occurs every time at least two of these wavepackets overlaps. The strength of the revival is 
determined by the number of wavepackets that overlap and these revival peaks occur at definite 
intervals. Using Eq. (|4ip we can write a general form for the revival times: 

k + -)t r , ( k+ — — ) t r , p = 1,2,..., N q and k = 1, 2, oo. (42) 



p J \ p 

Eq. (|4ip also highlights that when there are an even number of qubits the field returns to its 
initial state \a) after a time of t r . The time till the wavepackets return to their original position is 
determined by the slowest moving wavepackets, which are |e ±2i7ri / ir a) in the even qubit case and 
^ e ±mt/t r a ^ j n fo e odd qubit case. When N q is odd the wavepackets return to their start position 
after a time 2t r and there will always be no stationary wavepacket at the start position. 

A sketch of the Q function of the field for three qubits is shown in Fig. [9] for lar ge v alues of 



16| we can 



n and we once again show the dipole moments as arrows. Using \D k ) in reference 
compute a general term for the equation of the dipole moments in the large n approximation and 
is shown below. 

~~&k = (D k \ J N q \D k ) = (D k \~a Nq \D k ) 

^ lkl ( ±m ( e+ mii) i+si j e+ ?m)j) ( 43) 



where a N q = ®% q i + 3 + ®N q ^ ^ s * ne P aun vector for N q qubits and 

N q 

'* = $>f^. (44) 



~ x,y,. 
°N q 

i=l 



We notice again the pattern of the arrows changes over time and that the wavepackets come to 
their original starting position after a time 2t r . 

A. Basin of Attraction 

Using the above experience with one and two qubits, we investigate whether a basin of attrac- 
tion exists for N q > 2. Using the expansion pioneered by Chumakov et al. [is| and developed 
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FIG. 9: (color online) Phase space sketches of the Q function at eight different times when the qubits are 
initial started in the state \ggg). The atomic dipole states are represented as arrows, (a) the time t = 0. 
The system is started in a coherent state which is shown by a circle of uncertainty in phase space, (b) a 
time a little after t = 0. (c) the time t = (d) the time t = ^. (e) the time t=^f. (f) the time t = 
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independently by Meunier et al 



we have indeed found a basin of attraction for N„ qubits. As 



D_l_i(^ L )\ = l^tatt), and using the observation made in Eq. (jlU|) . we can therefore conclude that 
D,Nq_{t)) will go to a direct product of \^i att ) at t = jf^. 



By setting the values of (3 p (t) = 0, for p 



wavefunction and the qubits will go to 
a basin of attraction for N q qubits. 



®N q 



,<&N„ 



att 1 



■ — 2 then only 



31 2 



(45) 
(46) 
(47) 

(t) ) are in the 



2 2 

at i = Using this information we have found 



\V'. 



N q /2 

E 

m=-N q /2 



A{N q ,a,rn)e- i ^- m ^ e JWq~\ 



N q , m) 



N, 



No 



m 



A(N q ,a,m) 



«i 2 it(^ 



m is even 



m I is odd 



(48) 



As before <j> is an arbitrary phase and 6 is the initial phase of the radiation field with < \a\ < 
-j===. The basin of attraction is in the symmetric subspace and \N q , m) is the symmetrised state 
with N e excited qubits, N q qubits in the ground state and m 



N e -N a 



qubit case, I N q = 3, m 



For example, for the three 
4^ (\eeg) + \ege) + \gee)) as N e = 2 and N g = 1 etc. We note that 



although the state space of the qubits increases with increasing N q , the basin of attraction is still 
parameterised by the single complex value a. The significance of this will become apparent in Sec. 

We emphasise that when the qubits start in the basin of attraction then there are only two 
wavepackets in phase space of the photon field for all time. The initial conditions have restricted 
the system to behave like the one qubit case as is highlighted by considering the wavefunction at 

+ — 



-IOL, 



(49) 



2 x 2 

In the one qubit case the qubit system also goes to a second 'attractor' state \ipi att ) at the time 



^ and so 



,®N q 



att 1 







N q att 



. Qubits started in the basin of attraction, Eq. 



8|), will go to these 'attractor' states and hence the entropy will go to zero. 
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No. of qubits 


No. of Blobs 


Revival Times 


Possible 'Attractor' Times 


1 


2 








t r 

2 


2 


3 






tr , 
2 ' Lr 


tr 3tr 
4 ' 4 


3 


4 






3 ' 2 ' 3 ' r 


tr tr 5t r 

6 ' 2 ' 6 


4 


5 




tr 
4 


t r i r 2t r 3t r i 
3 ' 2 ' 3 ' 4 ' r 


tr 3r 5t r 7t r 

8 ' 8 ' 8 ' 8 


5 


6 


5 ' 


t r 

4 ' 


tr tr 2t r 3t r 4t r X 

3 ' 2 ' 3 ' 4 ' 5 ' r 


tr 3tr tr 7t r 9tr 
10 ' 10 ' 2 ' 10 ' 10 



TABLE I: A summary of where the different revival peaks occur for differing numbers of qubits. Also the 



times arc shown for when 



D ± ^ (t)) 



7V„ aft 



Note that only revival and 'attractor' times in the interval 



< t < t r are shown. Further times can be obtained by adding integer multiples of t r . 

We can construct a series for when these 'attractor' times will occur for N q qubits, shown in 
table [J with the revival times from Eq. (|42p . The general equation for these times is: 

(2p - 1) 



k + 



2N n 



p = 1,2,..., N„ and k = 1, 2, oo. 



(50) 



In Fig. [10] we show the 3 and 4 qubit evolution starting in a state in the basin of attraction (Eq. 
IS]) ). In general we do see there are dips in the entropy that go close to zero and the probability 
of being in the 'attractor' state goes to close to one. These features get more pronounced as n 
is increased. As for the two-qubit case, when this probability goes to zero, the system is in the 
second 'attractor' state, 



V'jv, att)> w hi cn is orthogonal to the first. 



B. Collapse and Revival of 'Schrodinger Cat' States 



Surprisingly, the basin of attraction has a simple form when written in terms of spin coherent 
states, which are analogous to coherent states of a photon field. As described in AppendixE] such 
states are given by: 

N q /2 i 



\*,N q ) 



1 



(1 + \ z f)N q /2 



E 



m=-N q /2 



N„\ (Nq , 

« z^^\N q ,m) 

(^-m)!(^ + m)! 



1 



T;(\e)+z\g)) 



(51) 
(52) 



(l + \ Z \ 2 )N q /2 

where \N q ,m) is defined in Sec. IIV Al A spin coherent state is a product state and this is easily 
seen when it is written in Major ana representation [40I in Eq. (|52p . Below we rewrite Eq. f|48j) in 
the form of a superposition of two spin coherent states, and show the form of the wavefunction at 
the time t = when the qubits are in the basin of attraction: 



*i\r 9 (0)) = |ViV 8 )|a> = K2^- 2 



1 



2 N " 



z = e- i9 ,N q 
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FIG. 10: (color online) Time evolution of the system when the initial qubit state is in the basin 
of attraction for n = 50 and 9 — 0. (Top) Three qubits: The initial state of the qubits is 



\N q = 3, m= |) + \/3\N q = 3,m = — ^))- (Bottom) Four qubits: The initial state of the qubits is 
ys(\N q — 4, m = 1) + \N q = 4, m= — 1}). As before (a) the entropy of the qubits. (b) the probability 



V2 

of all the qubits being in the state \g). (c) the probability of being in the qubit 'attractor' state 



N„ att 



+ 



Na-1 



2 N i 




(53) 



From Eq. (i52j) it is clear that \z = e l0 ,N q ) is orthogonal to \z 



-if) 



,N q ) and \^ Nq ) 



is a 



collection of spin 'Schrodinger cat' states. We also write the state of the field at the time of the 
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first dip in the entropy for an arbitrary number of qubits, t 



tr . 
2N„ ■ 



N a 



(&)) = K att ) [V¥^(a- ] [^-\aA^\ia) 



Nn-l 



2»i 



|,2 e -^n/2\_ ia \ 



(54) 



Eq. (|53p and Eq. (I54p are very similar and show how a Schrodinger cat state, and the infor- 
mation about the initial conditions of the qubits, moves from the qubits to the field. At the time 
t = the qubits are in a spin 'Schrodinger cat' state and the variable a, which fixes the initial 
state of the qubits, determines the ratio of \z = e~ iS , N q } to \z = —e~ ie , N q }. At the time t = 
the field is now in a 'Schrodinger cat' state and the information about the initial state of the qubit 
encoded by a is now stored in the field state. A 'Schrodinger cat' state, and the information about 
the initial qubit state has moved from the qubit system to the field system 34]. 



To illustrate this migration we use a function similar to the field Q function, Q q , using the 
completeness relation in appendix \A\ 



Q(a,t) 



(a\p f (t)\a) 



7T 



Qq(z,t) 



(z\p q (t)\ z ) 



(55) 



(i + M 2 ) 2 

Q q (z,t) corresponding to \^N q ) in Eq. (f53l) is depicted in Fig [Til There are indeed two separate 
peaks that correspond to two macroscopically different spin coherent states. 

The 'attractor' states for N q qubits in Eq. (|48p are also spin coherent states (see appendix |A|) . 
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(56) 



In Fig[T2]we plot a spin Q q function diagram for the 'attractor' state. 

So far we have shown, for large n, that the qubit states in the basin of attraction are spin 



Schrodinger cat states when a , ^ . At a later time, t = the qubits are in the spin 

coherent 'attractor' state and the field is in a Schrodinger cat state. Then the qubits return once 
again to a spin Schrodinger cat state at the time t = jj- . This is shown in Table [TTJ So we conclude 
that the 'Schrodinger cat' state migrates from the qubits to the field and back again. 
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FIG. 11: (color online) A plot of the spin Q q function of the qubits at the time t = for a state in the basin 
of attraction. The value of a = and we see a perfect Schrodinger cat state of the qubits. The number of 
qubits for the plot is N q = 40 with 6 = 0, but a similar picture can be produced for any N q > 3. The peaks 
merge together for N q < 3. 




FIG. 12: (color online) A plot of the spin Q q function of the qubits at the time t = t r /80, the 'attractor' 
time when N q = 40. The initial value of a = and we see a spin coherent state centered around the value 
z = ie lB with 6 = 0. 
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t = 


t = t r /2N q 


t = t r /N q 


Field 
Qubit 


Coherent State 
Spin Schrodinger cat state 


Schrodinger cat state 
Spin Coherent State 


Coherent State 
Spin Schrodinger cat state 



TABLE II: A table of the states of the field and qubits at different times. The behavior repeats in the large 
n approximation. 

C. Collapse and Revival of Entanglement 



Whilst the basin of attraction for the N q qubit case is a very small part of the Hilbert space 
of all initial states it shows many different degrees of entanglement, just like the two qubit case. 
Some values of a and 4> give the qubit state to be in a product state and other values give some 
degree of entanglement between the qubits. When considering more than two qubits there is no 
universally accepted measure of entanglement. Indeed, for N„ > 4 and multilevel systems there 



exist infinitely many inequivalent kinds of entanglement 



411 ] . GHZ type entanglement exists for 



all numbers of qubits and, up to local unitaries, a maximally entangled GHZ state has the form: 



\GHZ) = ± (\bf N « + 



(57) 



where {b\b±) = 0. In the previous section it was stated that the basin of attraction can be rewritten 
as a spin 'Schrodinger cat' state which is a state that is a superposition of two macroscopically 
distinguishable states. However expanding states in the basin of attraction using Eq. (|52p we see 
that it has a GHZ form: 
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(58) 



although it is only a maximal GHZ state when a = and a = —^==. Excluding the case when 
a = 2 nI-i i an other values of a give an entangled state of GHZ form. So, as with the previous case, 
there are only two values of a that give a product state when in the basin of attraction, meaning 
all other states contain some entanglement. We next consider the case of three qubits in more 
detail and verify by numerical calculations that the basin of attraction for three qubits is indeed 
only made up of GHZ type entanglement. 

For three qubits there are two different types of entanglement, GHZ 



421 and W 



4|. \GHZ) 



-j= (|eee) + \ggg)) can be regarded as a maximally entangled state, but if one of the qubits is 
traced out the remaining state has absolutely no entanglement present. The entanglement of 
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\W) = (\eeg) + \ ege) + |<?ee}) is maximally robust and is not so fragile to particle loss (i.e. 
there is still entanglement present after tracing out any individual qubit). It is not possible to 
transform \GHZ) into \W) even with a small probability of success and vise versa, and thus the 
two states represent fundamentally different types of entanglement. In particular, the GHZ state 
represents a truly 3-body entanglement. Similarly, any iVq-qubit GHZ state contains specifically 
./Vg-body entanglement. 



We have used a pure state GHZ entanglement measure introduced by Coffman et al. 



4|: 



TABC = C 2 



A(BC) 



AC ~ °AB 



(59) 



where tabc is the amount of GHZ entanglement, C^ BC ^ = 2y/detpA and C\ c is the entanglement 
after qubit B has been traced out. The results for different values of a for the three qubit version 
of Eq. (|48p are shown in Fig. [T3J This measurement only works for three qubit states that are 
pure, so we only perform the calculation on the initial state. 




FIG. 13: (color online) The value of the tangle for the states in the basin of attraction for different values 
of a. We notice that the plot shows only two points where the tangle is zero, just like the two qubit case. 



This plot shows the same shape as that of Fig. [6l but values of a are now limited to a smaller 
range. On the basis of Figs. [6land [T3"l and using Eq. ()58[) . we conjecture that for more qubits the 
basin of attraction will always contain maximally entangled states that involve all the qubits and 
two unentangled states, but the size of the basin will get smaller as the number of qubits increases. 

For the large n approximation it was shown in Sec. IIIIBI that if a state is started in the basin of 
attraction r(0) = r(%). This can also be extended to the three qubit case and ^^^(0) = tabc( gO- 
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Therefore if the initial qubit state is in the basin of attraction and tabc 7^ then the whole 
entanglement will return to the qubit system at t = |r. 

To see if 'Collapse and Revival' of entanglement also occurs in the exact solution for three qubits 
we start the qubits in an initial state, {1^3), that is a maximally entangled GHZ state (tabc = 1)- 
We then calculated the probability Pi n it{t) = (V , 3|Pg(£)|V'3)) which is the probability the qubits are 
in the initial state. If the probability collapses and revives over time then we can conclude that the 
three qubit system also exhibits 'Collapse and Revival' of entanglement. The results are shown in 
Fig. [Til The measure stated in Eq. [59] can not be used here as the qubits are mostly in a mixed 
state and Eq. [59] is for a pure state. For this reason we have chosen to look at the probability of 
being in the initial state. 

1 

0.8 
0.6 
0.4 
0.2 
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t/t r 

FIG. 14: (color online) 'Collapse and Revival of Entanglement' for three qubits for the initial state, 
\ (\[?>\m = 1/2) + \m = —3/2)), which is inside the basin of attraction, with n — 50 and 9 — 0. (a) 
the entropy of the qubits. (b) the probability of all the qubits being in the state \g). (c) the probability of 
being in the initial state, Pi n it(t)- (d) the value of tangle when one qubit has been traced out. In fact this 
pink line remains zero for all time, so there is no W state entanglement present. 

The qubits are started in a maximally entangled GHZ state, (tabc = 1)> that lies in the basin 
of attraction and are then allowed to evolve. At a later time the entropy goes close to zero and 
the qubits are close to the product form shown in Eq. (|47p . This is a completely unentangled 
state and as the entropy is also zero in the large n approximation, there is no entanglement in the 
entire system. After the 'attractor' time the entropy increases and is no longer close to zero. We 
follow the green curve in Fig. [141 and note that Pi n u(t) increases to a value near 0.8 and therefore 
we propose there is some three qubit mixed state entanglement present in the qubit system. In 
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the large n approximation this value does go to unity, so the system is in the initial state, up to 
a phase. In this approximation we will therefore necessarily get a full revival of three-body GHZ 
entanglement. 

As further confirmation that the entanglement is of 3-body form, we have also calculated the 
entanglement present between just two of the qubits after tracing out one of the qubits. It was 
found that the tangle between any two qubits remained zero throughout the evolution. Therefore 
we can conclude the entanglement present is just GHZ entanglement and there is no W component 
present (as a W component would lead to non-zero 2-qubit entanglement). 

We have demonstrated the revival of entanglement for the numerically evaluated exact solution 
for two and three qubits, and shown it to be true for all values of N q in the large n approximation. 
Of particular note is that the type of entanglement we observe is of the GHZ form and is thus a true 
iVq-qubit entanglement that cannot be reduced to a simpler form. Remarkably, this implies that 
the quantum information encoded in an A^-qubit GHZ entangled state — for any value of N q — can 
be encoded into the state of a single field mode. 



V. CONCLUSION 



We have studied the 'Collapse and Revival of Entanglement' in the two, three and man y q ubit 
JCM. We have shown that the idea of an 'attractor' state discovered by Gea-Banacloche [13j for 
one qubit does indeed exist in the multi qubit cases described by the Tavis-Cummings model. 
Unlike the one qubit case, the qubits have to be within a basin of attraction in order for the qubit 
system to disentangle itself from the field at certain times. When the qubits are started in a state 
in the basin of attraction then the system behaves in a similar way to the one qubit case. 

The basin of attraction contains mostly entangled states and there are some interesting entan- 
glement properties associated with the whole system when the qubits are started in this basin of 
attraction. The system evolves and there is then entanglement between the field and the qubit 
system, but when the system goes to the 'attractor' state all entanglement has completely left 
the system. The field and the qubits are both in pure states with the qubit state being in the 
completely unentangled 'attractor' state. At a later time it has been found that the entanglement 
does return to the system and then the process repeats. We represented explicit examples of en- 
tanglement present between the qubits and also between the qubit system and the field. We have 
called this phenomenon 'Collapse and Revival of Entanglement'. For initial qubit states outside 
the basin of attraction, sudden death and rebirth 3tJ of entanglement can occur. 
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APPENDIX A: SPIN COHERENT STATES 

There exist spin states analogous to the coherent state of the harmonic oscillator. For the one 
dimensional harmonic oscillator the coherent states are functions of a variable a which runs over 
the entire complex plane: 



\a) = le Qat |0) (Al) 
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»< - |0) (A2) 
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where the normalisation factor M = e' a ' / 2 . 

These states from a complete set, in the sense that j^j 

If °° 

- / \a)(a\d 2 a = ^ \n)(n\ = 1 (A4) 
77 J n=0 
where the right hand side is the unit matrix, or identity operator. 

To find the analogue for spin half particles we consider a single particle of spin S = N q /2. We 

define the ground state |0) as the state such that S z = S\0) where S z = Y^=o a i ls the operator 

of the z component of spin. Then the operator S + = sx + tSV creates spin deviations. In a simi 



way to the field coherent state the spin coherent state for many spin half particles is 



43, |46|,|47|: 



ar 



| Z ,JV,> = -^e'- S *\0) (A5) 

?1=0 

1 X 2 / m 



>Sit- m )\N q ,m) (A7) 



where the normalisation factor M = - — , \. N /0 , m = — n — is the spin of the state and z is 
a complex number that determines the center of the coherent state. A spin coherent state is a 
product state and therefore has no entanglement. For example if z = 1, then the average value of 
the spin coherent state is \N q , m = 0), if z = then the average value of the spin coherent state is 
N q ,m = ^jp ) and if z = 00 then the average value of the spin coherent state is N q ,m = — 



<q, — 2 

The states \z,N q ) form a complete set and the completeness relation is: 



N q + 1 f \z){z\ 
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In our work we present the 'attractor' states for many qubits in Eq. [56] as spin coherent states. 
These states themselves are indeed spin coherent states. If we set z = ie ld in Eq. IA7I then we get 



att ) an d ^ we se t z = —ie t9 then we get 



^Nq att 



When looking at the coherent state of a harmonic oscillator, there exists a set of states called 
Schrodinger cat states that are the addition of two coherent states that are macroscopically dif- 
ferent, e.g. -^i\ a ) + \~ a ))- There also exist spin coherent states that are analogous to these 
Schrodinger cat states. They are also defined as the addition of two macroscopically different 
states, i.e. -^(\z, N q ) + \-z, N q }). 

In our work we have plotted Q function diagrams for the cavity field state using the coherent 
state. There also exists an analogous function for the qubits, defined in terms of spin coherent 
states in a similar way, using the completeness relation Eq. IA8I This spin Q function is given by 

QqM = ^BM± 2 . (A9) 



l + \z\ 

We have used this in our analysis to make interesting plots which demonstrate pictorially that 
the states in the basin of attraction for multi-qubit systems are indeed Schrodinger cat states of 
the qubits. 
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